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The strong beams of high-frequency gravitational waves (GW) emitted by cusps and kinks of 
cosmic strings are studied in detail. As a consequence of these beams, the stochastic ensemble 
of GW's generated by a cosmological network of oscillating loops is strongly non Gaussian, and 
includes occasional sharp bursts that stand above the "confusion" GW noise made of many smaller 
overlapping bursts. Even if only 10% of all string loops have cusps these bursts might be detectable 
by the planned GW detectors LIGO/ VIRGO and LISA for string tensions as small as Gfi ~ 10"^'^. 
In the implausible case where the average cusp number per loop oscillation is extremely small, the 
smaller bursts emitted by the ubiquitous kinks will be detectable by LISA for string tensions as small 
as G/i ^ 10~^^. We show that the strongly non Gaussian nature of the stochastic GW's generated 
by strings modifies the usual derivation of constraints on Gjj, from pulsar timing experiments. In 
particular the usually considered "rms GW background" is, when G/j, ^ lO"'', an overestimate of the 
more relevant confusion GW noise because it includes rare, intense bursts. The consideration of the 
confusion GW noise suggests that a Grand Unified Theory (GUT) value Gfj, ~ 10~^ is compatible 
with existing pulsar data, and that a modest improvement in pulsar timing accuracy could detect 
the confusion noise coming from a network of cuspy string loops down to G/i ~ 10~^^. The GW 
bursts discussed here might be accompanied by Gamma Ray Bursts. 



I. INTRODUCTION 

Cosmic strings are fascinating objects which give rise to a rich variety of physical and astrophysical phenomena [Q. 
These linear topological defects are predicted in a wide class of elementary particle models and could be formed at 
a symmetry breaking phase transition in the early universe. Here, we shall reexamine the emission of gravitational 
waves (GW) by cosmic strings. The fact that oscillating loops of string are efficient GW emitters was pointed out long 
ago H] . The spectrum of the stochastic background of GW's generated by a cosmological network of cosmic strings 
ranges over many decades of frequency and has been extensively discussed in the literature ■ Until recently, it 
has been tacitly assumed that the GW background generated by cosmic strings was nearly Gaussian. However, in a 
recent letter prompted by a suggestion of 1^, we have shown that the GW background from strings is strongly 
non Gaussian and includes sharp GW bursts (GWB) emanating from cosmic string cusps. In Ref. we mentioned 
that kink discontinuities on cosmic strings also give rise to non-Gaussian GWB's. The simultaneous consideration of 
GWB's emitted by cusps and by kinks is important because, though cusps tend, generically, to form on smooth strings 
a few times per oscillation period |jl^, they might be absent on "kinky" strings (i.e. continuous, but non differentiable 
strings). On the other hand, the study of the process of loop fragmentation suggests that kinks are ubiquitous on 
loops (as well as on long strings) . 

In this paper, we shall discuss in some detail the amplitude, frequency spectrum, and waveform of the GWB's 
emitted both by cusps and kinks on cosmic strings. We shall also estimate the rate of occurrence of isolated bursts, 
standing above the nearly Gaussian background made by the superposition of the (more frequent) overlapping bursts. 
As we shall see, these occasional sharp bursts might be detectable by the planned GW detectors LIGO/VIRGO and 
LISA for string tensions as small as G/i ~ 10""'^^, i.e. in a wide range of seven orders of magnitude below the usually 
considered GUT scale G/xgut ~ 10"^. 

II. EMISSION OF GRAVITATIONAL WAVES BY COSMIC STRINGS IN THE LOCAL WAVE ZONE 

A. Waveform in the local wave zone 

In this Section we discuss the amplitude of the GW emitted by an arbitrary stress-energy distribution T^'^{x^) as 
seen by an observer in the "local wave zone" of the source, i.e. at a distance r from the source which is much larger 
than the GW wavelength of interest, but much smaller than the Hubble radius. For this purpose, we can consider 
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that the spacetime geometry is asymptotically flat: gjj'j^'*' = ry^i/ + hfj,u{x), where /i^j^ <C 1 is the metric perturbation 
generated by the source. The subsequent effect of the propagation of the GW in a curved Friedmann-Lemaitre universe 
will be discussed in the next section. 

Let us first consider a general scalar (flat spacetime) wave equation 

ULp{x,t)^-'iTTS{x,t), (2.1) 

and let us decompose the time variation of the source S in either a Fourier integral S{x, t) = J {dio/2-K) e~*"* S{x, lo) 
or (if the source motion is periodic) a Fourier series S{x, t) — e~'""* S{x, a;„). We concentrate on one frequency uj 
(or LOn). The corresponding decomposition of the solution, (p{x, t) — e^*'^* ^{x^ w), leads to a Helmholtz equation 

{A + uj'^)(p{x,uj) = -^■nSix.Lj) , (2.2) 

whose retarded Green function ((A + w^) Guj{x, x') — — 4:TtS{x — x')) is well-known to be Guj{x, x') — exp(-|- i wla; — 
x'\)/\x — x'\ so that 

/ iij \x-x'\ 
d'x' ^——^ S{x\u;). (2.3) 

If the source is localized around the origin {x' = 0) we can, as usual, replace, in the local wave zone {uj\x\ ^ 1), 
|a; — a;'| by r — n • x' in the phase factor, and simply by r in the denominator. [Here r = \x\, and n = x/r.] Let us 
define k = ujn (so that fc'' = {k,uj) is the 4-frequency of the (ys-quanta emitted in the n direction) and the following 
spacetime Fourier transform of the source 

S'(fc^) = S{k, = y d^x' e-'^-^' S{x', lu) . (2.4) 

With this notation the field if in the local wave zone reads simply 

(p{x,uj):::i Sik"), (2.5) 

r 



r 



where denotes either an integral over lu (in the non-periodic case) or a discrete sum over ujn (in the periodic, or 
quasi-periodic, case). 

Let us now apply this general formula to the case of GW emission by any localized source. We consider the linearized 
metric perturbation generated by the source: gft^{x^) = rj^i, -\- hfj,^{x'^). The trace- reversed metric perturbation 
hf^i, = h^^ ~ 5 ^ ^tJ-f satisfies (in harmonic gauge d'^ = 0) the linearized Einstein equations 

□ V ^ -ISTrGT^,, (2.7) 

where T^^{x^) denotes the stress-energy tensor of the source. We can apply the previous formulas by replacing 
Lp — > /ip^, S +4:GT^i,. Let us introduce the "renormalized" (distance-independent) asymptotic waveform k^^, 
such that (in the local wave zone) 

-h,Ax, t) = ^"'^^'""'"^ +o(\] . (2.8) 



Note the dependence of on the retarded time t ~ r and the direction of emission n. With this notation we have 
the simple formula (valid for any, possibly relativistic, source, at the linearized approximation [p^) 

K^,{t-r,n) = 4G ^ e~''^^'~^^ T^,{k,Lj) , (2.9) 

where we r ecall that k = ujn. In the case of a periodic source with fundamental period Ti, the sum in the R.H.S. 
of Eq. (2^) is a (two-sided) series over all the harmonics iw^ — ±muji with m e N and oji = 2tt/Ti, and the 
spacetime Fourier component of T^j/ is given by the following spacetime integral 

T;.(fc') = T^.(fe,c.) = ^ j\t I d?xe'^^''^-'^^T^,{x,t). (2.10) 
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B. Waveform emitted by a string loop 



We model the string dynamics by the Nambu action which leads to the string energy-momentum tensor 

T^^'ix^) ^fijdr daiX" X" - X'" X'") S^^'^ {x^ - X^{t, a)) . (2.11) 

Here, /i denotes the string tension and X^(t, cr) (to be distinguished from the spacetime point x^) represents the 
string wo rldsh eet, parametrized by the conformal coordinates t and a [X = dr X , X' = da X]. Inserting Eq. (2.11) 
into Eq. ( 2.10 ) yields the following Fourier transform of the string stress-energy tensor 

T^..{k^) = ^ [ dTda{X,X,-X'Xl)e-^''-'' . (2.12) 



Here k ■ X = rjf,^ k'' X" = X' - uX°, the indices of Tt"" have been lowered with r;^^ = diag (-f 1, -hi, -Hi, -1), 
and J2i denotes a strip of the worldsheet contained between two (center-of-mass) time hyperplanes separated by the 
fundamental period (denoted above as Ti) 

ive 2 

where ( denotes the "invariant length" of the closed loop that we consider. It is defined as I = E/ ji where E is the 
loop energy in its center-of-mass frame. [Note that i differs from the instantaneous length / da of the loop which 
changes as the loop oscillates.] 

The Nambu string dynamics in conformal gauge (and on our local, nearly flat, spacetime domain) yields a two- 
dimensional wave equation {&l — d^) X^(t, ct) = constrained by the (Virasoro) conditions 

r,^, X^ X'' + T]^, X'^ X"' = 0, 7?,,, X^ X"' = 0. (2.14) 

It is convenient to introduce the null (worldsheet) coordinates 

a± = T±a, j^ = d±^l-{dr±da), (2.15) 
oa± 2 

and to decompose in left and right movers (note the factor i) 

X^ir,a)^^[Xtia+)+X^_ia^)]. (2.16) 

In terms of this decomposition, the Virasoro conditions read (9+ X^)"^ = = (9_ X^)"^. We can (and will) also choose 
a (center-of-mass) "time gauge" in which the worldsheet coordinate r coincidcf^ with the Lorentz time in the center- 
of-mass frame, i.e. X^{t, a) = r, so that X^ = cr+, = o'_, and X]. is (for a closed loop in the center-of-mass frame) 
a periodic function of a± of period £. In this time gauge, the Virasoro conditions yield (X+)^ — 1 — (X_)^ where 
the overdot denotes the derivative with respect to the corresponding (unique) variable ct+ or cr_ entering X_|-((t+) or 

Inserting (2.16) in (2.12) yields the following result for the Fourier transform of (to be inserted in the waveform 

(H) 



where X^X'^^ = ^ (X^ -I- X'^ X'^) denotes a symmetrization on the two indices fiv, where is a truncated 
cylinder on the worldsheet defined, say, by < r < = £/2, and < a < £, and where we recall that 



'^Note that the requirement that = r (without any proportionality factor) links the period P± in a± to the value of the 
string center-of-mass energy E, namely: P± = £, with E — fi£. In the fundamental-string literature the periodicity in a± is 
fixed to be, say, 27r (for a closed string), and one writes = t/2tv fi — Et jl-K ^ = £t/2'k. 
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A 47r 

Kn = {k.m,UJrn) = {ni UJg U, m IVg) = — TO(n, 1) 



(2.18) 



with m e Z — {0}, runs over the discrete set of the 4-frequencies of the GW emitted by a string of invariant length £ 
in the direction n. [In the following, we shall sometimes restrict m to the positive integers, it b eing understood that 
one must then add the complex conjugate quantity when computing the asymptotic waveform (2.9).] 

The result ( 2.17| ) can be further simplified by changing the variables o f inte gration from (r, a) to {<tj^ , cr_ ) . One 
must use dr da = -k da< da_ and take care of the limitation of the integral (2.17) to the truncated cylinder y^„. This is 



( 2.17 ) to the truncated cylinder J^t 
most easily done by rewriting (2.17) as /i£ times the average over the worldsheet ^ /J"* dr /Tg da/lj of the integrand 

X- exp —ik ■ X. Remembering that the period in a± is £, the averaging can be rewritten as J^da^/£ J^da-/£. 
If we postpone the symmetrization on the indices fj.v to the last stage of the calculation we can write 



" £ ^ 



where we introduce the following asymmetric double integral 



d<74 



+ 

Jo 



(2.19) 



(2.20) 



Using the complete factorization of the integrand of (2.2C) in the product of a function of cr+ and a function of a-, 
we can finally write 



(2.21) 



where 



da±Xte-i''--^± 



The final factorized (modulo the symmetrization) result 



(2.22) 



(2.23) 



0) follows from the easily 



will be very convenient for our subsequent study. The conservation of T^'^ (i.e. 
checked identity km^Jt = satisfied by the simple integrals ( 2.22 ). 

Note that Eq. (|^ gives a factorized expression for the Fourier transform of the GW amplitude ( |2.9| ). Such 
left-right factorized expressions are ch aract eristic of quantum amplitudes of closed (fundamental) string processes. 
The convenient factorized expression ( p. 23 ) was used in our letter Q for the calculation of the classical radiation 
amplitudes of cosmic closed strings in the Fourier domain. Previous calculations of GW amplitudes were performed 
in the time domai n |13yil[| , though the factors (2.22) appeared as building blocks in the calculation of the radiation 
power from loops ]14|-|16|. The publication of our work Bthen prompted other authors to recognize the convenience 
of left-right factorization in GW amplitude calculations |17| . 



C. Decay with frequency of the waveform: cusps, kinks and other singularities 



If we define 4>±{o'±) = ki- X±{a±) (where ki is the m = 1 v alue of (2.18) ), the hi gh-fre quency behaviour of T^^{k )^ 
and therefore of the Fourier transform of the waveform ( |2.9| ), is reduced, by Eq. ( 2.23 ), to studying the asymptotic 
behaviour, as m ^ oo, of the two simple integrals I±{m) = da± f±{<J±) with f±{(J±) — X!^{a±) e-5™**±('^±). As 
is well-known the asymptotic behaviour (as m oo) of l!^(rn) depends on essentially two things: (i) the regularity 
(i.e. the number of continuous derivatives) of the functions X!^{a±) and (/)+((T+), and (ii) the presence or absence of 
saddle points (stationary-phase points) in the phase 0±((T±) (i.e. of points (tJ where 0±(o'±) = 0)- If the hmctions 
X!^{a±) are smooth (C°°) and if <j)+{<^+) and 0_(ct_) have no saddle points the integrals I±{m) tend to zero faster 

than any negative power of m as rn ^ oo. A fortiori, the product T'^'^{k) oc I^(m) l'^{m) then tends to zero faster 
than any negative power of m. In such a case, the GW emission of one string loop would be well approximated by 
considering only a few of the lowest mode numbers m. 
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By contrast, in the present paper we focus on the case where (i) and (ii) are violated in such a way that T^^{k) has 
a rather slow, power-law decay as to — > oo. The two physically most relevant cases where this happens are near cusps 
or kinks. First, note that (say) the + phase 0-)_((t+) = -q^^, /cf X^((T+) has a saddle point 0+ = when fci • X+ = 0. 
Remembering that both and (because of the Virasoro condition) are null vectors, we see that saddle points 
occur each time fc^, and therefore k^^ = m k'^ , is parallel to X^^{a^). In the time gauge, where = {X±, 1) with 

— 1, Jl + ((t+) and Jc_(cr_) correspond to two separate curves, say C+ and C_, on the unit sphere [^ . The 
saddle points occur when the unit direction vector n of the emitted GW lies either on C+ or C_. If one has only 
one saddle point, say in the phase factor of /^(to), the integral /^(to) will have a slow decay as to — > oo. But if the 
other integral I^(m) has neither a saddle point, nor some lack of regularity in X^{a^), the integral If^{m) will decay 

exponentially fast with to, so that the product T'^'^{k) cx l'^\m) l'^ (m) will still decay exponentially fast. 

Therefore, the two generic cases where T'^'^{k) can have a slow, power-law decay are: (a) the case where the two 
curves C+ and C_ intersect, so that the fc'^ parallel to their intersection develops a double saddle point, or (b) the case 
where k'^ is parallel to a direction of C+ (or C_) and where the dual function X^{a^) (or, respectively, X^{a+)) has 
some type of discontinuity. The case (a) corresponds to a cusp, and happens generically for smooth (and in particular 
continuous) closed curves C± PQ I . The discontinuity in the case (b) can be of various type (say a mild discontinuity 
in some higher derivative of X^{(t±)). The most interesting case (leading to the slowest decay with to) is the case 
of a kink, where, say, X^(cr+) is continuous but X^(cr_) has one or several jump discontinuities. It is expected that 
kinks are ubiquitous on loops (and on long strings) . Note that the presence of kinks (which is expected because of the 
reconnections) means that the curves C± on the unit sphere are discontinuous. [Hence, too many kinks can prevent 
the two curves C± from intersecting, i.e. can prevent the presence of cusps 



D. Logarithmic Fourier transform of GWB waveforms 



For the time being, we wish to conclude from this discussion that, in the presen ce o f cusps or kinks, the discrete 
Fourier components of the asymptotic wave form K^i^(u;m,Ti) oc T^i,{km,iiJm)i Eq. ( |2.9D , will decay as to — > cx) in a 
slow, power-law manner along certain directions: a finite set of directions (corresponding to the intersections of C+ 
and C-) in the case of a cusp, or a one-dimensional, "fan-like", set of directions (corresponding to C+ (and/) or C_) 
in the case of kinks. If the observer at infinity happens to lie near one of those special directions, it will detect a 
stronger than usual GW amplitude: these are the gravitational wave bursts (GWB) that we study in this paper. 

We see that, by definition, the GWB's correspond to a large value of the harmonic number to, i.e. to a frequency 
fm = Wm/27r = m/Ti = 2m I i much larger tha n th e frequency of the fundamental mode of the string. For such 



high mode numbers to the discrete Fourier sum (2.9) can be approximated by a continuous Fourier integral (indeed. 



Acj — i^rn+x ^ i^rn = i^i = ^m/m ^ Wm). In Other words, on the time scales At of relevance for the detection of 



GWB's (/„! < At < T^) we can replace in Eq. (2.9) 



so that 

K^,{t~r,n)^2Gl j ^ e--(«-'-) T^^(fc, ^) . (2.25) 

To any continuous function, say n{t), of some (possibly retarded) time variable t we associate the following logarithmic 
continuous Fourier component k(/) (corresponding to an octave of frequency around the analyzing frequency /): 

K[f)^\f\n{f)^\f \ I dt e'-^f' K[t) . (2.26) 



[The advantage of this definition over the straightforward Fourier tran sform /?(/) is that k(/) has always the same 
physical dimension as K[t).] In terms of this definition, the result ( ^.25| ) leads to the following simple formula for the 
logarithmic Fourier transform of the GWB asymptotic waveform: 

ti^,{f,n)^2Gl\f\T^,{k,Lo). (2.27) 



When inserting the factorized form ( 2.23 ) this yields more explicitly 
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'if,n) = 2G^,\f\Ii^r_} 



(2.28) 



where the simple integrals were defined by Eq. 



Remember that these expressions give the asymptotic (distance-independent) waveform (2.8) in the local wave 
zone of the source, and that the frequency / still refers to the frequency measured in the local wave zone of the 
center-of-mass frame of the source. The problem of the cosmological propagation of will be discussed later. 



III. GRAVITATIONAL WAVE BURSTS EMITTED BY CUSPS AND KINKS 

A. Waveforms from cusps 

As recalled above, a cusp corresponds to an intersection of the two curves C+ and C_, i.e. to a point on the 
worldsheet where (in the time gauge) the two null vectors X^{a+) and X^(cr_) coincide. Let us denote 



r = in^^\ 1) = = X^iaf) 



Ac). 



(3.1) 



the common value of these two null vectors at the cusp X^^^ = X^((T^'' , ctI^-*). The (spacetime) direction of strongest 
emission from the cusp is precisely i.e. the GWB is centered around the 4- frequencies fc^ oc i.e. remembering 
Eq. ( ^l8| ), the space direction of strongest emission is n — n^'^'. Let us first study the Fourier transform of the 
waveform emitted precisely at the center of the GWB (i.e. n = n^^\ and fc^ — muj^i^^ ). We shall discuss below the 

(c) 

beam width around this direction. To simplify the writing we shift the origin of a± so that CTj_ = 0, and the origin 
of so that — 0. We can then write the following local Taylor expansions (truncated to the order which is 
crucial for our purpose) 



(3.2) 



X^{cT±) - + a± + ^ Xf^ a\ 



(3.3) 



where the successive derivatives (with x'^^ = dj. X±) appearing on the R.H.S. are all evaluated at the cusps (i.e. at 



0). Differentiating the Virasoro constraints X^. 



Therefore, at the cusp, one has i ■ X± = and £ ■ X. 



(3) 



yields the relations X± ■ X± ^ and X± ■ X^^' + X^ = 0. 
— {X±)'^. [From which one sees that Xt^ is a spacelike vector.] 



These rela tions yield the following simple result for the crucial quantities k ■ X± oc £ ■ X± entering the phase factor in 
Eqs. (^ItI) or (p2|) 



e,xf^{a±) = -- (X^)V^ 

[This shows a posteriori why it was c rucial to include the terms ©((t^) in the local Taylor expansion of Xf^{a±) 
Inserting these results in Eq. (2.22) leads to an expression of the form 



(3.4) 



da±{i^' + X^a± + ^ 



(3.5) 



As we shall see the intervals of cr+ and cr_ which contribute most are (because of the saddle point in the phases) very 
small for large m (Acr± oc |m|~^/^). It would then seem that the dominant term in l!^ is obtained by keeping only the 
leading term in the parenthesis, i.e. i^^, so that l!^ ~ a±£'^. However, this leading contribution does not correspond 
to a physical GW, but can be removed by a coordinate transformation. Indeed, as we are working in the Fourier 
domain (and with the asymptotic waveform), a linearized coordinate transformation has the following effect on k^i,: 



(3.6) 



Here, we are considering the case fc^ oc As k'*'' oc if we decompose = a± i^^ + b^., where bj_ denotes the 

subleading contribution from (3.5), both the leading- leading term a+ a_ i'^ , and the two leading-subleading terms 
a+£'^b'^ and b'^a-H" can be gauged away. [This explains why our final waveform below differs from that obtained 
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earlier in Ref. 



which did not notice that the leading terms were pure gauge.] Finally, the leading, physical 



waveform is given by keeping only I^I_ 



b^fb''_} with 



ao+e 



d a± a± exp y — muji X^ cr^ 



(3.7) 



[We used fc^ = mujgi^, where we recall that oji = ^ir/Ti ~ 47r/i' is the basic loop circular frequency, linked to the 
GW frequency by / = l o JI-k = moJi/i'K = Imj i with m e Z.] 

Most of the integral (^J) comes from a small interval in a± around zero. This allows us to neglect the limitation 
to a period [croi fo + ^] a-nd to formally extend the integration on from — oo to + cx3. It is convenient to introduce 
the scaled variables 



u± — N± a± ; N± 



1 
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This leads to the appearance of the following integral (the same for 71+ and u_) 



/ = 



(3.8) 



(3.9) 



Here, the sign ± denotes the sign of m, i.e. the sign of the frequency /. 
by an interval of order unity in u = u±, corresponding to Aa± ~ ^/N±. 
pure imaginary and to be 



27r 



3r I 



It is clear that the value of / is dominated 
The exact value of / is easily found to be 



(3.10) 



where F denotes Euler's gamma function. Note that the square of /, which enters the waveform, is real, negative and 
independent of the sign of m. Finally, if we define 



X^ 

\xJl 



47r(12)^ 



(3r(i)) 



2 ' 



(3.11) 



we find that the (logarithmic) Fourier transform of the asymptotic waveform reads (for positive or negative frequencies) 



C 



(27r 1/1)1/3 



(3.12) 



Here, we have introduced the arrival time of the center of the burst, ic, which was set to zero in the calculation 
above (by our convention X^^^ = 0). The fact that the two (generically independent) spacelike vectors A^. (which are 



orthogonal to i'^) are real means that the GW (3.12) is linearly polarized. 

To understand the meaning of the |/|~"'^^'^ dependence of the Fourier amplitude ( |3.12 ), we take the inverse Fourier 
transform (remembering the definition ( p. 26 )) which yields a time-domain waveform proportional to 



K(t) CX |t — • 



(3.13) 



Note that the fact that Eq. ( |3.13 ) tends to zero a,t t — tc do es no t mean that the GWB is best detected away from 
t = tc- The full waveform, in the time-domain, is the sum of (3.13) and of a slowly varying component (due to the low 



modes of the string). What is important, and distinguishes the GWB from the slowly varying component, is the fact 
that (3.13) is "spiky" , because of the appearance of the absolute value of i — tc- If one were to consider the curvature 
(tidal forces) associated to ( [3.13| ) it would be oc |i — <c|~^, exhibiting more clearly the spiky nature of the GWB. 

Actually, the sharp spike at t — tc exists only in the limit where the observer lies exactly, at some moment, along 
the special direction n*^^^ defined by the cusp velocity, i.e. when n = n^^^. Let us define 6 as being the angle between 
the direction of emission n and the "3-velocity" of the cusp n''^^ We shall now show that when ^ 9 <^ 1 the 
time-domain cusp waveform is approximately given by the 6 = waveform computed above, except in a time interval 
around tr of order 



t-tc 



0'Tf 



(3.14) 
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where the spike is smoothed. In the frequency domain this smoothing on time scales (3.14) corresponds to an 
exponential decay for frequencies 



l/l^^- (3-15) 

To study the effect of 7^ 0, let us introduce the four vector (5'' such that l'^ — fc'' + S'^ where = (n, 1). In the 
time gauge 5'^ = (n'"^) — n, 0) is spacelike and of squared norm — 2{1 ~ cos 6') ~ 9^. Therefore (5^ = 0{6). Going 



back to the expression (|3.5|), and remembering from Eq. (3.6) that one can gauge away any term in k^i, having a 



factor kf^ (X kfj_, we see that we should now split the parenthesis in Eq. (3^) as k^ + {S^ + X^^ cr± . . .) and decompose 
accordingly II^ — a±k'^ + bj. with 



o"oH 



b'i= dermis" + X^a± + ■■ ■)e-^"'"^^ . (3.16) 



By a gauge transformation we can, as above, discard the a± k^ term and replace I'^ by hj^. On the other hand, in the 
phase terms we have now (using k ■ £ = — | (fc — £)'^ = — ^ ~ — ^ 6^ and k ■ X± = —6 ■ X±) 



K X^{(T±) = k-ia± + ]^k-X±al + ^k-Xf al 



-i ^?V± - 1 <5 . X± 4 - 1 {X^f al , (3.17) 



instead of (B.4) 



If we rescale a± as in Eq. (|3.8|) and introduce 



e± = = 9 ' '.. ' (3.18) 

\X±\ \12\X±\) 

we see that the gauge-simplified value of (i.e. ( |3.16| )) is, after factorization of an overall factor ~ \X±\/N'^, of the 
form (when neglecting factors of order unity, and treating and 5^ as scalars) 

^±(£±)= /dw(e± +w)e''^±("^^±) , (3.19) 



where (t)±{u,£±) ~ + e± + £\u. Remembering that the integral ( ^.9| ) is dominated by what happens in an 
interval Au ~ 1, this schematic expression is sufficient for seeing that when ej. ^ 1 the numerical value of /±(e±) is 
well approximated by /±(0) = /, i.e. that K^y(/, n) ~ K,xjy(/, 't-^'^^)- To discuss what happens when, on the contrary 
e± > 1 one must study a little bit more carefully the behaviour of the phase 4)±{u, £±) ^ + £±v?' + £\ u. Let us go 



back to the unsealed expression (3.17) and differentiate it: 



-J-(k,X^{a^))^\5^ + {5,X>^)a± + \{X±fol. (3.20) 

The discriminant A = 6^ — 4ac of this trinomial in (t± is A = (5 • X^f' — S^{X±)^ — —S^iX^)^ sin^ /3± where /3± is the 
angle between the two space vectors 6 and X±. The important point is that (generically) A < which means that 
the trinomial ( 3.20| ) has no real roots, i.e. that 0±(cr±) has no saddle point when 6* 7^ 0. In fact, this absence of saddle 



point when 6* 7^ can be seen, as an exact result, from the fact that in the scalar product k^J_X^Jya±) both 4-vectors 
are null and future-oriented so that their product can vanish only if they ar e par allel, but we wanted to show how. 



within certain limits for the unwritten numerical coefficients, the toy integral (3.19) with '/i±(u,£±) ~ "u?" -\-£± u^+e^. u 
could qualitatively represent the exact result for all values of £± (both > 1 and <1). The absence of saddle point 
means that when £± gets significantly larger than one, /±(e±) will tend exponentially fast toward zero. As there are 
only numbers of order unity in the (unwritten) coefficients of I±{£±) we conclude (as usual for such estimates) that: 
(i) when e± < 1, /±(e±) can be estimated by /±(0) — I (though this estimate is numerically accurate only if e± ^ 1), 

while (ii) when e± <^ 1, /±(£±) starts decaying exponentially fast with £±. Consequently, the waveform 
will also interpolate, as 9 increases, between essentially «'*''(/, n^^^) and an exponentially small result. 
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As neglecting factors of I-k might be detrimental to our subsequent estimatesg we tried to be a little bit more precise 
about these orders of magnitude. First, let us note that the facts that (in the notation used here) the period in (j± 
of XjJyG±) is P± — i and that X± are unit vectors imply that the generic order of magnitude of \X±\ (if the string 
is not too wiggly) is 



(because X.± — c„exp(27rz7i(T±/i?)). Using the estimate (|3.21), using a± 



(3.21) 

^ and writing that the divide 



between small 9's and large 9's is obtained when the third term on the R.H.S. of ( 3.17 ) is equal to the first leads to 
6»dividc = (4/(V3^|/|))^/^ = (2.31/£|/|)i/3. Approximating 2.31 ~ 2 leads to the simple resuh 



idc 



i2/\f\i)'^' = i\m) 



-1/3 



(3.22) 



where Ti = i/2 is the basic period of the string motion. This corresponds to the inequality (3.15) quoted above. 
When passing from the Fourie r domain to the time domain, the exponential decay in the domain (3.15) (now consider 
for a fixed 0, instead of (3.22) which considered / as fixed and let vary) means that the waveform becomes smooth 
on time scales At 6^ Ti near the center of the GWB, as announced in Eq. (3.14) above. 



As we a re di scus sing " 27r-accurate" estimates, let us conclude this subsection by mentio ning t hat when inserting the 
estimate ( 3.21 ) in ( ^.12 ) there appears the coefficient (12)3 /^/(27r^) (with /im given by ( 3.10 )) which is numerically 
= 0.8507, i.e. close enough to 1 to be neglected. Finally, a good estimate of the amplitude of the asymptotic waveform 
(when one is not interested in polarization effects) is simply 



K™'^P(/, n) 



_Gjxi_ 

(l/K)^ 



6(61 

divide (/) 



-i(n-n(=))) 



(3.23) 



where 8(2;) is the step function (1 if a; > 0; if a; < 0). One should remember that the result ( |3.23 ) has been derived 
by assuming that |/| ^ = 2 |m| ^ 1. As the asymptotic GW amplitude generated by a string at low frequencies 
I/I t — 0(1) is 0(0)1 1), we see that, amp litude- wize, a cusp GWB is a small correction to a low-frequency background. 
But what is essential in the result (3.23) is the very slow decay with the mode number, oc |r7i|~^/^. 



B. Waveforms from kinks 



GW emission by kinks has been studied by Garfinkle and Vachaspati |ll||. However, like in the case of cusps, the 
leading term that they studied turns out to be pure gauge. This will be clear from the different, Fourier-domain, 
treatment that we give now, which is a simple generalization of the method discussed for cusps in the previous 
subsection. 



As discussed above, kink emission corresponds, in the original expression Eq. (3.17), to the case where, say, the 
phase 0+(cr+) = fci • X+ has a saddle point (or is close to a saddle point), and where Xt{cr-) has a discontinuity (at 
some CT_ — a^^'^). Though the discussion is somewhat different than for the cusp case, we shall be brief as the method 
of attack is a variant of the one we discussed in great detail above. The saddle point requirement for implies that 

fc^ must be nearly aligned with some null vector X^{a^^'^). [As said above, the set of all exactly aligned null vectors, 
i.e. the set of all the central null geodesies within the beam emitted by a moving kink, is a one-dimensional, fan-like, 
structure defined by fixing (T_ = a'^^'^ , and letting (t+ run over its entire period.] The most convenient starting point 

is again the factorized form of T^''(fc^) cx I^l'^ , where we recall, for convenience, the form of the simple integrals 



J (Jn 



(3.24) 



where we introduced a gauge parameter ^± whose value can be (and will be) chosen to simplify l!^. The integral 



I'^ is treated as in subsection II A above (using some ^-|- 7^ 0) with the same results (including the effect of 9 



^We mention this because the "surprisingly large" value of the parameter F < 
a loop can essentially be attributed to a factor in F. 



50 entering the total rate of GW energy loss of 
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COS ^ n ■ n^'^^ 7^ 0). In particular, we recall that (after gauging away some terms) the value of in the aligned case 
is 

l1~XlJda+ C7+ 6*"'^'^+'"+ , (3.25) 

which scales with m (as m — > ±oo) like ±|m|~^/'^ — \m\^^^^m^^^^ (where ± is the sign of m). 

On the other hand, calls for a new treatment. In fact, if we assume that X'^{a-) jumps from X^((t1""^ — 0) — rii 
(a null vector) to X'^ (ctI"''^ + 0) = rij (another null vector), we get the leading estimate of the integral by replacing 
X^(cr_) by n^' (cr_ - crl'"'^) for cr_ < (tI'"'^ and by (cr_ - crl'"'=) for cr_ > ctI'^'^ (we set X^(crl'"'=) = 0). [We are here 
following standard results on "edge effects" in oscillatory integrals, see, e.g., |l^.] This yields (in the m ^ oo limit, 
and with ^_ = 0) 



~ 2i ~ . (3.26) 

' K • ni k ■ n2 



The essential feature differentiati ng th e result ( 3.26 ) from the normal "cusp" result ( |3.25 ) is its scaling with ra as 



^±00. The kink contribution ( |3.26 ) decays as m ^ = ±|to| ^, i.e. faster (by |m| ^/'^) than the ±.\m\ decay of 

( 3.25 ) . When considering the waveform , we can finally contrast, in order-of-magnitude, 

the previous "cusp" result, Jj^^^^ = /^'^"^p /^'^"'^p to the new "kink" one J^l'^y. = /^'^"''p /^kink_ Therefore the ratio 



.kinky^cusp _ jkink/jcusp esscutially given by the ratio /kink/jcusp^ j ^ ^j^^ ^.^^j^ between (|3.26|) and the usual 



result (3.25) (with H > — ). If the discontinuity in is of order |^-| '^-^ 1, the ratio /kmk^jcusp giniply given 

(independentl y of t he sign of m) by the power |m|~^/'^ ~ characterizing the faster decay of the (simple) 

kink integral (|3.26| ) versus its cuspy analog. This simple reasoning allows us to i mmed iately translate our previous 
cusp results into their kink analogs. The — i power of |/| in the cusp waveform ( 3.12| ) becomes replaced by a — | 
power: 



^k^„k(/)«l/r^e2"/*=4"5-^ (3-27) 



with B'^ proportional to the {6 = Q limit of) the vector Eq. (i.2i). The time-domain waveform becomes 

K'''"'^(t)cx (3.28) 



and still corresponds to a formally infinite spike in tidal GW forces as t ^ tc- Finally, the simplified estimate ( 3.23| ) 
translates into 

^''"'(^' ~ e(0divide(,/) - cos-i(n . n(*^))) , (3.29) 

where n^*^^ is the direction closest to n within the "fan" radiated by the moving kink. 

Let us note that our general, Fourier-domain approach can easily deal with weaker types of GW emitting worldsheet 
singularities. For instance, if we consider a weaker kink where X^(cr_) and X'^{a-) are continuous, but where 
X^((T_) is discontinuous, the decay of ( ^.26 ) will be replaced by a to ~^ de cay as m — * 00. This faster decay will 



correspondingly increase (by one) the (inverse) power of |/| appearing in ( 3.27 ) and ( 3.29] ) 



Let us also note that Ref. [g/7| has recently studied the waveforms emitted by piecewise-linear loops, i.e. the case 
where both X'^ and are piecewise constant, with discontinuities at some kinks. In this very special case both 
and I'^ are given by a finite sum (over the number of kinks) of terms of the form ( 3.26| ) (corresponding to one 



kink), in which one must reinsert the kink phase factor exp(— z/c • X_/2). The scaling with m corresponding to this 
case is Jl^'^^^^ — /^'^'"'^ J^^mk ^ ^-2 _ |r7i|-2^ This leads to a waveform oc \'fn\I^_^y.^^]^l''J ^^^^y. oc |m|~^. The 

time-domain version of this k''"''^'^(/) cx |/|^^ is, near each kink, K''"°^''(i) cx |< — tc|^^- The corresponding time- 
domain curvature vanishes everywhere, except at the discrete set of kink arrival times where it has a delta-function 
singularity. We thus recover the finding of that the time-domain waveform of such piecewise-linear loops is 
a piecewise-linear function of retarded time. Our analysis shows, however, that such special piecewise-linear loops 
are bad models of the waveforms emitted by generic string loops. Indeed, even if a string network contains only 
a small fraction (say a few percent) of loops with cusps, this small fraction will dominate (see below) the crucial 
high-frequency tail of GW emission (because k™''p(/) oc \f\-^'^). Even in the a priori implausible case where the 
fraction of cuspy loops is negligibly small, the high-frequency tail of GW emission will be dominated by generic 
kink waveforms (esc -f+^cusp-^-^kink)' with K^™^{f) oc |/|~^^'^). The faster decay of the special piecewise-linear loops, 
^linear °^ l^^l'^+'kink^- kink ^ I/I '^^'^ disqualifies their use as models of GW emission by a network of strings. 
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IV. PROPAGATION OF A GRAVITATIONAL WAVE BURST IN A COSMOLOGICAL SPACETIME 



In the previous Sections we discussed the emission of a GWB in the local wavezone of the source, i.e. at distances 
large compared to the wavelength but small compared to the cosmological scale. T he G WB amplitude was then 
characterized by its (distance-independent) asymptotic amplitude k^i/, entering Eq. (|2.8| ). We need now to study 
the subsequent effect of the propagation of h^i, in a cosmological spacetime. It is well-known that if we consider a 
perturbation, g^i, = -t- h^^, away from an arbitrarily curved background spacetime gj^^{x^), the trace-reversed 

perturbation /i^i, = h^i, — ^ g^^^ g^°'^ hap satisfies, in the gauge hap = 0, and away from the source, the propagation 
equation 

5^ Kp + 2 h^^ - 2 i?^^ /i^)^ = , (4.1) 

where denotes the covariant derivative defined by the background metric. We consider the case where our GWs 
have wavelengths much smaller than the scale of variation of the background metric. [This is certainly the case for 
what concerns the cosmological background. We shall not consider here the special situations that arise when the 
GW meets, during its propagation, some local bump in the curvature, of scale comparable to its wavelength.] In 



such a case, we can: (i) neglect the curvature terms in Eq. (4.1), and (ii) treat the leading propagation equation 



9b' ~ in the WKB approximation: 



/ia;3 = Real [Aca/se*^/^] , (4.2) 

where the polarization tensor is normalized by g"^^ g^ Cap e^u — 1. 

As usual the WKB approximation {e 0) yields, if we introduce the wave vector = S/e (with k'^ = g^ k^)^ 



9b ^fj. ^i' — , 



(eikonal equation) 



v^(A2 kn 



0, 



A{2k^' A + fc^ A) = 0. 



(4.3a) 
(4.3b) 
(4.3c) 
(4.3d) 



For our present purpose, the most important results are Eqs. (4.3c) and (4.3d). Eq. (4.3c) says simply, in words, 
that the transverse (see Eq. ( 4.3b| )) polarization tensor Cap of the GW is parallely propagated along the n ull ge odesies 
(Eq. (4.3a)) describing, in the geometrical optics limit, the GW propagation. Most important is Eq. (4.3d) which 
gives the law of decrease of the G W am plitude A along the null ray. 

If we write down the condition ( 4.3cj ) for the case of a spatially flat Friedniann-Leniaitre universe 



ds^ = -dt^ + a^{t){dr^ dn^) = 0^(77) [-drf + d?^ + dn^ 



(4.4) 



and for a "retarded" solution of the eikonal equation (4.3a) of the form S = F{ri — r), (where we choose the center- 



of-mass of the source as center of the polar coordinate system) we find that r a A remains conserved during the 
propagation, i.e. that the GW amplitude decreases as 



A = 



0(77) r 



(4.5) 



77— r— const. 



In the local wave zone a{r])rc^ o- iVcrn ) ^ ~ f (where the subscript "em" refers t o th e emission event) is the physical 
radius r which appeared in Eq. ( ^.8[ ), so that the constant k on the RHS of (4.5) measures the amplitude of the 
asymptotic GW tensor amplitude k^j/, after having factorized the normalized polarization tensor e^i/. Finally, the 
time-domain GW amplitude arriving on Earth can be written as 



(^rcc) 



(4.6) 



where tree denotes the proper time at reception, arcc = a(^roc), ??roc = / dt/a{t), and where "pp" means that the 
tensor must be parallely propagated, between the emission and the reception, along the null geodesic followed by 
the GW. As the latter null geodesic is described by dt/a{t)~r — 0, we have the usual redshifting of time intervals 
between emission and reception, dtrec/flrec = dtem/acm, which corresponds, in the Fourier domain, to /rocOrec = 
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/cm = (1 + /roc , 1 + Z = a^oc/ O-cm ■ (4.7) 

The logarithmic^ Fourier transform of the GW amphtude at reception can be written in terms of the fogarithmic 
Fourier transform of the asymptotic GW amphtude at emission, Kfii,{fcm), as 

KPP((l + z)f) 

^ ■ (4-8) 
ao r 

Here, and henceforth, / = /rcc denotes the observed frequency, ao = arcc denotes the present cosmological scale factor, 
and z the cosmological redshift introduced in Eq. ([1.7|). It remains to express the "amplitude distance" aor (which 
is (1 + z)^^ times the luminosity distance) in terms of the redshift z. We use for this the relation valid in a spatially 
flat, matter-dominated {flmo = 1) universe: 

aof= 3.0(1-^), (4.9) 

where to — 2/{3Hq) is the present age of the universe. [In the numerical estimates below, we use Hq ~ 
65 kms~^ Mpc~^ which corresponds to to — 1-0 x 10^" yr ~ 10^^'^ s.] Though this relation gets modified in the 



earlier radiation-dominated era, it will be sufficient for our purpose to use Eq. (4.9) for all values of z, because agr 
tends anyway to the finite limit Sto = 2/Hq as z gets large. 



In the following, we shall work with order-of-m agnit ude estimates. We simplify the "am plitude distance" (4.9) to 
aoT to z/{l + z), and use our simple estimates ( 3.23| ) (for the cusp waveform), and ( |3.29 ) (for the kink waveform). 



Finally, we have (in terms of the observed frequency / ~ /rcc, henceforth considered as being positive) 

((l + z)/^)V3 toz ' 

and 

Note that the low frequency part ((l + z) / ^ T^^ ^ £^^, i.e. low mode numbers |r7i| ~ 1) of the GW amplitude would 
be of order h^^ ~ Gjil/aor ^ Gfj,£{l + z)/{toz). Compared to this non-burst, "full" signal, we have the simple 
orders of magnitude : /i'="''p(/) ~ 6l,„(/) /i^F ^nd ft'^'"'=(/) ~ 0l^{f)h^^, where 0^(1) ~ {{I + z) f e)-^/^ - |to|-i/3 
embodies the basi c pow er-law dependence on the mode number m when |m| 3> 1. It is crucial to keep in mind that 
the "cusp" result (P^ holds only if, for a given observed frequency /, the angle 9 between the direction of emission 



n and the 3-velocity n*^^-* of the cusp satisfies 

'<0„,= ((l + z)/^/2)-i/3. (4.12) 



some 



Similarly, the "kink" result ( |4.1l| ) holds only if the sma llest angle between the direction of emission n and 

kink velocity 3- vector n'^*^) satisfies the same relation (4.12). Note that the domain of validity of the cusp result ( |4l0| ) 



is, for each loop period, a (small) cone, of half opening 0,n, around n^''\ while the domain of validity of the kink 



result (4.11) is a 0,„-thickening of the "fan" of directions drawn by the continuous time evolution of the kink velocity 
vector TV^. 

V. GRAVITATIONAL WAVE BURSTS FROM A COSMOLOGICAL NETWORK OF STRING LOOPS 

A. Simplified description of a string network 

Having derived the GW amplitudes emitted by individual cusps and kinks on some loop situated at cosmological 
distances, we need now to sum the contributions coming from a cosmological network of string loops. For this, we 
shall use a simplified description of such a string network. Indeed, though much work has been done to understand the 



Note that the definition (2.26) ensures that a constant redshift affects the argument of k(/) but not its amplitude. 
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evolution of such networks (see references in |^), there remain many uncertainties about some of the crucial detailed 
features of this evolution (notably the exact value of the parameter a introduced below, and the average number of 
cusps per loop). In fact, our work provides a new motivation for reinvestiagting such questions and getting better 
answers. Anyway, in the present exploratory investigation we shall content ourselves with using a very simple ("one 
scale") description of a string network. Let us recall that, at any cosmic time t, a horizon-size volume contains a few 
long strings stretching across the volume, and a large number of small closed loops. The typical length and number 
density of loops formed at time t are approximately given by 

£r^at , ni{t) a^'^ t^^ . (5.1) 



As we said above, the exact value of the (crucial) dimensionless parameter a in (O) is not known. We shall assume, 
following 1^, that a is determined by the gravitational backreaction, so that 

a-TG^, with r- 50. (5.2) 

The coefficient F is defined as that entering the total rate of energy loss by gravitational radiation d£/dt — TG fi^. 
[Note that fundamental string theory suggests that string loops of small size loose energy not only as gravitons, but 
also as dilatons, which increases the effective value of F [Q.] For a loop of invariant length £ (and oscillation period 
Te = e/2) the lifetime is e/TG^j,^ t. 

In the following, we shall express all the cosmological dependence in terms of the redshift z, rather than the cosmic 
time t. Let 

Zcq ~ 2.4 X 10"* nrnO K ~ lO^'^ (5.3) 

denote the redshift of equal matter and radiation densities. For z < Zoq, i.e. during matter domination, we have 
a{t)la^^{tlt^fl^^{l + z)-\i.e. 

t ^ta{l + zY^/"^ (matter era). (5.4) 

On the other hand, for z > z^q (radiation era) we have (1 + z)^^ = a{t)/aQ — (acq/ao)(t/tcq)^^^ so that 

t^toil + z,qy^^l + z)-\ (5.5) 

For our subsequent estimates, we found convenient to define smooth functions of z which interpolate between the 
diffe rent functional dependences of z in the matter era, and the radiation era. For instance, in view of Eqs. (5.4) and 
(5.5) we define the smooth function 

Mz) = (1 + z)--'/^ (1 + z/zeq)-'/' , (5.6) 

in terms of which, 

t C:i to ipi{z) . (5.7) 



Then, from Eq. (5^), the typical length of a loop formed (and decayed) around the redshift z is 

i^atoMz), (5.8) 

while their number density is 

a-H^^ipJ^iz). (5.9) 



B. Gravitational wave bursts from cusps 



In this subsection we concentrate on cusp GWB's. Inserting Eq. ( p.8\) into Eq. (4.10) yields a GW amplitude from 
cusps at redshift z of the form 

/i™^P(/, z) ^ G^^a'/'{fto)-'/' M^) 6(1 - &„.(«, /, z)) , (5.10) 

where we defined the interpolating function 
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^n{z) = z-\l + (1 + z/zeq)-'/' , (5.11) 

and where the 0-function factor {Q{x) denoting as above the step function: G(x) = for a; < 0, 0(a;) = 1 for 
a; > 0) serves the purpose of cutting off the burst signals that would formally correspond to 0„j > 1. Indeed, the entire 
derivation of the burst signal in Section |l| was done under the assumption ^ 1, corresponding to high values of 
the mode number |m| ^ 9:^. The low mode numbers m = 0(1) do not correspond to bursts, and the string does not 
emit modes with |m| < 1. We mentioned above that h^^^^{f) ~ dm{f)h}"^ where h^^ ~ G fj,£/aor is the amp litude 
of the low frequency signal generated by the low mode numbers. Therefore, formally the cusp signal ( 5.10| ) gives 
an approximate representation of the string GW amplitude which is valid for all Omif) > 1, down to and including 
the (formal) limit 9^ = 1- The explicit expr essi on for 9m{ct, f , z) is obtained from combining Eq. ( 4.12| ) (where we 
henceforth neglect the factor 2^/^) with Eq. (5.8), and reads 

0™(a, /, z) = {af toil + z) Mz))-'^'' = (« / to)-'^^ (1 + z)^/'' (1 + z/z,^)^/'' . (5.12) 

Let us now turn to the problem of estimating the rate of occurrence of GWB's from cusps. As recalled above, for 
smooth loops cusps are generic, and tend to be formed a few times during each oscillation period pO[ |. Reconnection, 
and its associated kink formation, can, however, diminish the average number of cusps |1^. We find it, however, 
very plausible that a significant fraction of the loops will exhibit cusps. To quantify this, we introduce a parameter c 
defined as the (ensemble) average number of cusps per oscillation period of a loop. 

We start by estimating the rate of GWB's originating at cusps in the redshift interval dz, and observed around the 
frequency /, as 

dN- iC(l + z)"V(z)dy(z). (5.13) 



Here, the first factor is the beaming fraction within the cone of maximal angle 9mif, z), Eq. ( 5.12 ); the second factor 
comes from the link d^obs = (1 + 2) dt between the observed time iobs (entering the occurrence rate on the L.H.S.) 
and the cosmic time t of emission; the quantity 

i,(t)^^^^2ca-H-^ (5.14) 

is the number of cusp events per unit spacetime volume (in which enters the average number c of cusps per loop 
period Tg = £/2 ^ at/2); and, finally, dV{z) denotes the proper spatial volume between redshifts z and z + dz. In the 
matter era, 

dV = 54 7rtjj [(1 + z)^/2 - if (1 + z)-"/^ dz , (5.15) 

while in the radiation era 

dV = 72 Trig (1 + Zcq)^/^ (1 + z)-^ dz . (5.16) 

It is convenient to work with the logarithmic density N{f,z) = dN/d In z. Using the relations given above, we 
write it in terms of a new interpolating function of z: 

N{f, z) ~ 10^ c tg-i a-8/3 (/ to)-2/3 ^^^(^) ^ (5 17) 

where the numerical factor 10^ approximates an exact numerical factor which is 54 7r/4 when z < 1, 54 tt when 
1 < z < Zeq, and 72 tt when z > Zgq, and where we defined 

ip„iz) = z3(l + z)-V6 (1 + ^/^^q)ii/6 . (5.18) 

The observationally most relevant question is: what is the typical amplitude of cusp-generated bursts h^'^^^{f) that 
we can expect to detect at some given occurrence rate N, say, one per year? As the function (/?„(z) always increases 
with z like a power-law (with an index which depends on the considered range of redshift), the value of TV is dominated 
by the largest redshift, say Zm, contibuting to N: 

N{f, z) d In z ~ N{f, z„0 . (5.19) 
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The looked for estimate h}^'^^^{f) is therefore obtained by: (i) solving Eq. (5.19) for z™, or, equivalently, solving 



Eq. (5.17) for z, and (ii) substituting the result z — z,n{N, /) in Eq. (5.10). The final answer has a different functional 
form depending on the magnitude of the quantity 



j/(7V,/)^10-2(iV/c)toa«/=^(/to)'/^ 



(5.20) 



Indeed, if y < 1 the dominant redshift will be Zm{y) < 1; while, if 1 < ?/ < ?/eq = zlq^^, 1 < Zm{y) < Zcq, and if 
y > j/cq, Zmiy) > Zcq. More precisel y, th e solution of Eq. (5.17) for z can be written as the following (interpolating) 
function of the combination y, Eq. ( 5.2C| ): 



z^iy) = y'/Hl + yy^''{l + y/y.q)-'/'' 



(5.21) 



above. 



where yeq 

We can again introduce a suitable interpolating function g{y) to represent the final result as an explicit function of 

y- 

\f)^Gfia'^'{fto)-'^'g[y{NJ)]e{l 



^cusp f 



(5.22) 



where 



5(2/)^y-i/3(l + 2/)-"/33(l + y/yeq)^/" 



(5.23) 



where Q{x) denotes as above the step function, and where ^,n(a, N, /) denotes the function of a, N and / obtained by 
substituting z Zm{y{N, /)) (defined by Eqs. ( 5.2C| ) and ( ^.2l[ )) into Eq. (5.12) above. In fact, this 8-function cutoff 
will be needed only when we consider very low frequencies / and very small values of a. For instance, if / ^ 1/ (7yr) 
and N/c ^ 1/yr, O mja, N, /) would become larger than one only for a < 10~^. 

The prediction ( ^.22 ) for the amplitude of the GWB's generated at cusps of cosmic strings is one of the central 
results of this work. Before proceeding to analyzing the detectability of these bursts, let us discuss the GWB's 
generated at kinks. 



C. Gravitational wave bursts from kinks 



We recall that, from Eqs. (4.11), (4.12), the two differences between kinks and cusps are: (i) the kink GW amplitude 
is smaller than the cusp one by a factor \m\-^^^ ^ Om ^ ((1 + z) f i)'^/^ (i.e. h^"'^ - 6^^ h^^ instead of /i'^""? ~ 
(^mh^^), and (ii) the kink amplitude is emitted (per period) in a thickened fan of directions of solid angle ~ 9m, 
instead of a cone of solid angle ~ This second fact is in favour of the kink signal, but we shall see that it does 
not suffice to com pensa te the bad news that the kink signal is parametrically smaller than the cusp one. 

Using formula (5.12), we can easily derive the kink analogues of the cusp results derived above. First, we find, 
instead of ( 5.1C ) , 



^kink 



(/) 



(a, /, z) /i™^P(/) ~ Gfia'/^ if to)--'^-' ip^^' (z) 9(1 - 0„(a, /, z)) , 



2/3 ,Ak) 



with the kink analog of the cusp interpolating function (5.11): 



{k) 



(z) = Z-l(l + z)-l/6(l + z/Zeq)-'/^ 



(5.24) 



(5.25) 



The rate of GWB's originating from kinks in the redshift interval dz, and observed around the frequency / is obtained 
by dividing Eq. ( [sl^ ) by em{aj,z), Eq. ( |5T2| ). This yields, instead of Eqs. ( p7| ), (|t|). 



^^kinks/^ In z ~ 102 a-^'^{fto)-^/^ (z) 



where 



«(z) = z3(l + z)-4A^(l + z/ze,)'V3 



(5.26) 



(5.27) 



The parameter k in Eq. ( ^.26 ) is the kink analogue of the parameter c in (5.17), i.e. the average number of kinks on 
a loop. Now, the expectation is that A: > 1. In the following we shall simply assume A: ~ 1, though one must keep in 
mind that k might be significantly larger than 1. 
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Like in our discussion of cusps, we are interested in estimating the GW amphtude of k ink b ursts that one can 
expect to detect at a given recurrence rate N. As before, this is obtained by first solving Eq. ( 5.26| ) for z, which yields 



(5.28) 



where ycq^ 



"cq"^, and where the quantity y''^' is the following function of N and /: 



y« (TV, /) EE 10-\N/k) to (ji^)i/3 ^ 



(5.29) 



We can finally write 



'kink 

'■N 



(/) ~ G M a'/' if to)-'/' [4'^ b^'^ (N, /)]] 6(1 - el^^ («, TV, /)) , 



(5.30) 



where oln* (a, N, f) denotes the result of substituting z 
(^.30) in terms of an interpolating function g'^'^^ (jj'''^^), as in Eq. (5.23) 



(y W{N /)) into Eq. (|5.12|) . We could also have written 



D. Functional behaviour of /i'^^^p and h^'"^, and first comparison with planned GW detectors 



It is easily checked that both ft-™^'^(/), Eq. (5.22), and (/), Eq. (5.30), are monotonically decreasing functions 



of both N and /. Note also that /i™**^ depends on the average number of cusps c only through the combination 
N/c, while h^'^^ depends on the average number of kinks k only through N/k. The dependence on N and / (as 
well as c and k) can be described by (approximate) power laws, with an index which depends on the relevant range 
of dominant redshifts. Let us focus on the functional dependences of /i™''^, which will turn out to be the physically 
most relevant quantity. [It is easy to derive the analogous results for /i^"*^ by using the formulas given above.] As 
N increases (or as c decreases), h'^^^^ decreases first like N^^^^ (or c^/"^) in the range Zm < 1, then like N~^^^^ 
(or c®/^^) when 1 < Zm < Zeq, and finally like N~^^^^ (or c'/^^) when z„j > Zcq- For the frequency dependence 
of h'^^^P, the corresponding power-law indices are successively: —5/9, —9/11 and —7/11. [These slopes come from 
combining the basic f~^/^ dependence of the spectrum of each cusp-burst with the indirect dependence on / of the 



dominant redshift Zfyi{a, N , f), Eqs. (5.21), (5.2C).] By contrast with these monotonic behaviours, when using our 
assumed link G fi ^ a/50 between the string tension and the parameter a, one finds that the index of the power-law 
dependence of h'^^^'^ upon a takes successively the values: -1-7/9, —3/11 and +5/11. The appearance of the negative 
index —3/11 means that in a certain intermediate range of values of a (corresponding to 1 < z^(a, N , f) < Zoq or 

1 < y{a,N,f) < ycq — Zcq^^) the GWB amplitude (paradoxically) increases as one decreases a, i.e. G fi. [A decrease 
of a leads to a smaller radiation power from individual loops at a given redshift, but at the same time it also leads 
to a higher density of loops and thus to a higher likelihood for an observer to see some of the loops at a small angle 
with respect to cusp direction. The overall effect is determined by the interplay of these two factors.] 

In Fig. |l| we plot (as solid lines) the logarithm of the GW burst amplitude, log]^o(^*'"^*); ^ function of log]^o(Q^) 
for: (i) cusps with c — 1 (upper curve), (ii) cusps with c = 0.1 (middle curve), and (iii) kinks with k = 1 (lower curve). 
Fig. uses the fiducial value N ~ 1 yr~^, and gives the value of h™^'^ or h^"^^ for a frequency f = fc = 150 Hz. As 
is discussed in the next Section, this central frequency is the optimal one for the detection of a /"^^'^-spectrum burst 
by LIGO. We indicate on the same plot (as horizontal dashed lines) the (one-sigma) noise levels /i"™*''' of LIGO 1 (the 
initial detector), and LIGO 2 (its planned advanced configuration). The VIRGO detector has essentially the same 
noise level as LIGO 1 for the GW bursts considered here. We defer to the next Section the precise definition of these 
noise levels, as well as the meaning of the extra short-dashed line in the lower right corner of Fig. |l|. 

From Fig. |l| we see that the discovery potential of ground-based GW interferometric detectors is richer than hitherto 
envisaged, as it could detect (if c ^ 1) cosmic strings in the range a > 10"-^°, i.e. G iJ,> 10~^'^ (which corresponds to 
string symmetry breaking scales > 10"'^'^ GeV). Even if c ~ 0.1, i.e if cusps are present only on 10% of the loops in the 
network, which we deem quite plausible, (advanced) ground-based GW interferometric detectors might detect GW 
bursts from cusps in a wide range of values of a. Let us also note that the value of a suggested by the (superconducting- 
) cosmic-string Gamma Ray Burst (GRB) model of Ref. |9|, namely a ~ 10^^, nearly corresponds in Fig. 0to a local 
maximum of the GW cusp amplitude. [This local maximum corresponds to 1. The local minimum on its right 

corresponds to z ~ Zeq.] In view of the crudeness of our estimates, it is quite possible that LIGO 1/VIRGO might 
be sensitive enough to detect these GW bursts. Indeed, if one searches for GW bursts which are (nearly) coincident 
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FIG. 1. Gravitational wave amplitude of bursts emitted by cosmic string cusps (upper curves) and kinks (lower curve) in the 
LIGO/VIRGO frequency band, as a function of the parameter a = 50G n (in a base-10 log-log plot). The upper curve assumes 
that the average number of cusps per loop oscillation is c = 1. The middle curve assumes c = 0.1. The lower curve gives the 
kink signal (assuming only one kink per loop). The horizontal dashed lines indicate the one sigma noise levels (after optimal 
filtering) of LIGO 1 (initial detector) and LIGO 2 (advanced configuration). The short-dashed line indicates the "confusion" 
amplitude noise of the stochastic GW background. 



with (someQ) GRB the needed threshold for a convincing coincident detection is much closer to unity than in a bhnd 
search. Let us finally note that Fig. ^ indicates that (except if k happens to be parametrically large) the kink bursts 
are too weak to provide an interesting source for LIGO/VIRGO. 

In Fig. ||, we do the same plot as Fig. |l| (still with N — 1 yr~^), but with a central frequency f = fc ~ 3.9 x 10^'^ Hz 
optimized for a detection by the planned space borne GW detector LISA. The meaning of the various curves is 
the same as in Fig. |l|. The main differences with the previous plot are: (i) the signal strength, and the SNR, are 
typically much higher for LISA than for LIGO, so that LISA could be sensitive to even smaller values of a (down to 
a ~ lO^^^'®), (ii) LISA is very sensitive even to rare cusp events (c = 0.1 or even smaller), (iii) LISA is, contrary to 
LIGO, sensitive to the kink bursts (which are believed to be ubiquitous), and (iv) though the GW burst signals still 
stand out well above the cusp-confusion background (discussed in the next Section) , the latter is now higher than the 
(broad-band) detector noise in a wide range of values of a. LISA is clearly a very sensitive probe of cosmic strings. 
We note again that a search in coincidence with GRB's would ease detection. 



VI. DETECTION ISSUES, CONFUSION NOISE, PULSAR TIMING EXPERIMENTS 

A. Signal to noise considerations 

Let us first complete the explanation of Figs. ^ and |^ by discussing the choice of the central frequencies and the 
detector noise levels indicated there. 

We recall that the optimal squared signal to noise ratio (SNR) for the detection of an incoming GW by correlation 
with a suitable bank of matched filters is given by 

Here h{f) is the Fourier transform of the (best) template (assumed to match the signal), Sn{f) is the (two-sided) 
noise spectral density, and, as above, we introduce the logarithmic Fourier quantities h{f) = |/| h(f), and hn{f) = 



■^The local maximum of the in Fig. |l| corresponds to a redshift Zm. ~ 1. By contrast, in the model of ^ the (300 

times more numerous) GRB's come from a larger volume, up to redshifts ~ 4. 
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FIG. 2. Gravitational wave amplitude of bursts emitted by cosmic string cusps (upper curves) and kinks (lower curve) in 
the LISA frequency band, as a function of the parameter a — 50G n (in a base-10 log-log plot). The meaning of the three solid 
curves is as in Fig. The short-dashed slanted curve indicates the confusion noise. The lower long-dashed line indicates the 
one sigma noise level (after optimal filtering) of LISA. 



^nif)- For cusp bursts (that we focus on) the optimal bank of fihers (when neglecting the fine structure around 
the center of the cusp) is 



h{f) = A\f\-'/^ 



(6.2) 



and depends, besides the overaU amphtude factor, on only one parameter: the arrival time tc- We take the foUowing 
model of the LIGO 1 (two sided) noise curve (see, e.g., ||l^) (for / above the seismic cut off fs ~ 40 Hz) 



Snif) 



■ Sq 



2 + 2 



(6.3) 



and /o 200 Hz. [Th e for m (|6.3[) is not real ly up to date, but it is sufhcient for our 



By inserting Eqs. (3.2) and (3.3) into Eq. (6.1) we get an explicit integral proportional 



Here, 5*0 = 1.47 x lO-^^Hz" 
present orientation estimate. 

to f{dx/x){l/s{x)) where x = ///o and s{x) — x^^{2 + 2x' + x~^). The minimum of the function s{x) is located 
at Xm = 0.7483, which corresponds to fm = Xm fo — 149.67 Hz. Therefore LIGO 1 is optimally sensitive, for such 
signals, to the frequencies / ~ fm ~ 150 Hz. [This value would also be approximately appropriate for kink signals, 
and also for the LIGO 2 noise curve.] Choosing fc — 150 Hz as fiducial frequency, and reexpressing the full SNR (6.1) 
(including its overall amplitude factor cx \A\'^) in terms not of A but of h{fc), one finds (after computing the integral) 
that 



with an "effective" noise level 



Loff 



MM 



1.7 X 10 



-22 



(6.4) 



(6.5) 



The effective noise level ( |6.5[ ) (which corresponds to a SNR = 1 for a matched filter detection^) is what is called the 
one Sigma noise level" of LIGO 1 in Fig. |l|. For LIGO 2 (advanced configuration) we estimated from noise curves, 
available on the LIGO web site, that, near 150 Hz, the noise amplitude hn{fc) is a factor ~ 13.5 smaller than for 
LIGO 1. This defines the lower dashed curve in Fig. ||. 

We did a similar analysis for LISA. We used as (effective) noise curve the sum (with a factor 1/2 included to take 
care of our using a two-sided spectral density) 



^The effective noise level (6.5), corresponding to the horizontal line in Fig. ^ (with a similar line in Fig. should not be used 
to estimate the SNR for the detection of a stochastic background, which is optimized by a different filtering technique. 
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si°'{f)^l[sr''{f) + s'fr'{f)], (6.6) 



where S'jf''*'' is a recent estimate of the instrumental contribution to the noise |2l| 



Sr''if) ^ 2.13 X 10-^^ 



f ' ^' 

J a 

7 



) , (6.7) 



where fa — 2.76 x 10 ■^Hz, fp — 9.55 x 10 ■^Hz, and where 5'™"^(/) is the "binary confusion noise", as estimated in 
Ref. Again the optimal frequency is fixed by considering the minimum of ^^(/)/|/i(/)P oc /^^^ 5'^°*(/), which 

occurs at = ]^q-2.4ii3 _ 3 379 x 10~3 JJ2. And again the "one sigma" effective noise level (dashed horizontal 



line in Fig. ||) is defined by (|j), with the result: hf = 1.815 x 10^22 

Let us also briefly mention the problem of "thresholds", i.e. the minimum value of the SNR, say po, needed to 
distinguish, with enough confidence, a real signal from a statistical fluctuation of the instrumental noise. Assuming, 
for simplicity, Gaussian instrumental noise, the probability that the template-filtered detector output exceed a certain 
level Po of SNR is given by the complementary error function 

p{p > po) - ' dpe-i"'' ^ 0) ' (po - pI) e-^^o . (6.8) 

We are interested in the situation where two independent detectors (either two ground based interferometers, or the 
two, partly independent, subinterferometers of LISA) find a coincidence, after having made N observations during 
the year. In our case, each template contains only the arrival time as essential free parameter. Therefore, the number 
of observations, assuming a blind search, in one year is TV ~ (1 yr)/T = (10^'^ s)/t where r ~ is the characteristic 
time between two decorrelated successive observations. Assuming the same noise lev el in each detector, the looked for 
threshold po can be defined by equalling the product of two probabilities equal to ( |6.8| ) (one for each detector), i.e. 
the square of (|6.8|) , to 1/N. When r ~ 10"^ s (as appropriate to LIGO) this yields po — 4.3, while when r ~ 3 x 10^ s 
(LISA) this yields po ~ 3.0. Note, however, that the model of ^ suggests that GWB's may be associated to Gamma 
Ray Bursts. The threshold for a search in near coincidence with GRB's is somewhat lower, because of the smaller 
number of trial observations N. 



B. Confusion noise due to gravitational waves from a string network 



The realization that the stochastic ensemble of GW's generated by a cosmological network of oscillating loops is 
strongly non Gaussian, and includes occasional sharp bursts, raises the following issues, of crucial importance for the 
detection strategies: (i) Can one split this stochastic ensemble of GW's into a (strongly non Gaussian, but plausibly 
nearly Poissonian) "burst" part (best detected by a matched filter approach) , and a (nearly Gaussian) "background" 
(best detected by the usual strategies discussed for Gaussian stochastic backgrounds)?, and (ii) What is the relation 
of this split to previous estimates of the "stochastic" string-generated background of GW's l^-jj , and how does it 
affect the interpretation of the famous pulsar timing constraint |2^? 

Our proposal is, for each detector with characteristic (optimal) detection frequency /c, and for each GW amplitude 
level, to define the borderline between occasional, individual sharp bursts, and a nearly Gaussian background by 
counting the average number of bursts of given amplitude which arrive within a characteristic time tc = fc^- In. other 
words, we define a nearly Gaussian background by considering the confusion noise generated by the overlap of more 
than one (and generally many) bursts which arrive within a time smaller than the considered characteristic inverse 
frequency. A technical way of justifying this (physically intuitive) consideration is the following. 

Let us write, in the time domain, our stochastic ensemble of GW's from a string network as 

h{t)=Y,Kt~in.Zr,,Pn). (6.9) 
n 

Here, the (somewhat symbolic) sum runs over waveforms arriving at time t„ and emitted from a string loop at redshift 
z„, w ith other string parameters (length, orientation,. . .) being denoted p„. The (logarithmic) Fourier transform of 



( p. 9]) reads 

h{f) = I/I h{f) - e'"^*" '^(/' ^"'^^") ■ (6.10) 
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Let us now recall that the spectral noise density of a stochastic ensemble of signals is defined by {h*{f)h{f')) = 
d{f — f')Sh{f) where ( ) denotes an ensemble average, a tilde the usual Fourier transform and a star, complex 
conjugation. Defining the (total) root mean square (rms) GW amplitude /irms(/) (with the same dimension as h{t), 
i.e. dimensionless) by h'^^^if) = \f \ Sh{f), the above definition of Sh{f) becomes, in terms of the more convenient 
logarithmic Fourier quantities 



(/i*(/)M/'))-l/l^(/-/')/^Ls(/)- 

When we take the (formal) limit /' /, the delta function S{f — /') 
infinite) total time interval T = J dt 

{h*{f)h{f)) = \f\Thl,M)- 



(6.11) 

J dtexp{2TTi{f — f')t) becomes a (formally 

(6.12) 



Let us now compute the quantity {h*{f)h{f)) by squaring the expression ( |6.10|) . Before taking the ensemble 
average, we get a double sum over n and n' involving phase factors e^'^'-''^*" in addition to the other phase factors 
hidden in the dependence on the other parameters. If we assume (as usual) that such phase factors are random and 
average to zero (except when n = n') we get 



{h*if)hif)) =Y,\Hf,^n,Pn 



3.13) 



Within our simplified approach to the string network, we assume that the GW amplitudes differ only by their redshift 
of emission z„. The sum over z„ (within some octave around z) then counts the number of signals coming from dz/z 
during the total time T. In terms of our previously introduced differential rate of occurrence N{ f, z) this yields simply 



{h*{f)h{f))=T 



dz 



N{J,z) h^if.z) 



Identifying this result with ( 6.12 ) wc finally get 



dz 



htmsif)^ I -niLz)h'{f,z). 



where we introduced the shorthand notation 



N{f,z) _ 1 dN 



(6.14) 



(6.15) 



(6.16) 



This derivation has achieved two aims: (1) it gives us an explicit expression, ( |6.15| ), (computable in terms of 
quantities that we calculated above) for the usually considered rms GW background generated by a string network, 
and (2) it shows (by comparison to the usual rms value of a sum of n independent random variables with the 
same variance) that the quantity n{f,z), Eq. (6.16), gives, in a technically precise sense, the (effective) number, 
within an octave of frequency around /, of random GW bursts generated at redshift z, and therefore of amplitude 
h{f,z), which contribute to h^^^. The latter result leads us to split the ensemble of GW signals in two sets: (i) 
the set of rare, non overlapping bursts such that n(f,z) < 1, and (ii) the set of frequent, overlapping bursts, such 
that n(f,z) > 1. The rare, non overlapping bursts contribute to h^^g{f) only if one considers integration times 
T ^ (iV(/, z))^^ = (n(/, z) l/l)"'^ > |/|~^- Therefore, if we are interested in detection issues involving a detector 
with a certain characteristic bandwidth ~ /c, and a corresponding integration time Tc ^ fc^i the bursts such 
that n{fc, z) < 1 should be considered as randomly occurring sepa rate burst events, and the magnitude of these 
occasional events should no t be compared to the full fi^^^ifc) of Eq. (3.15) but only to the "confusion" noise defined 
by restricting the integral ( 6.15 ) to the overlapping events, n(/c,z) > 1. [By the central limit theorem, applicable 
when n{fc, z) ^ 1, this confusion noise can be considered as being nearly Gaussian.] Therefore we define 



confusion 



(/) 



dz 



n{f,z)h\f,z)Q{n{f,z)~l) 



(6.17) 



We conclude that the relevant background that individual cusp or kink bursts should exceed to be detectable by 
LIGO (central frequency fc = 150 Hz) or LISA (/, = 3.88x lO'^Hz) is not /iL.(.fc), Eq. ( |6l5| ), but only /i?onfusion(/c), 
Eq. ( |6.17D . The short dashed lines in Fig. ^and Fig. |^ plot precisely the quantity (6.17), for the c = 1 cusp background, 
i.e. for h{f,z) given by ( |5.10| ), and N{f,z) given by Eq. ( |5.17| ) (with c = 1). This shows that in the LIGO or LISA 
bandwidths the individual bursts occurring once per year stand out clearly above the relevant confusion noise. 
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C. Rare bursts, confusion noise and pulsar timing experiments 



Our finding that the stochastic ensemble of string-generated GW's is not Gaussian, but can be viewed as the 
superposition of occasional bursts on top of a nearly Gaussian "confusion" background leads us to reexamine the 
pulsar timing experiments |23| and their use as constraints on the string tension G jjL. Let us take as characteristic 
frequency of the pulsar timing experiments the frequency /p*^"^ = 1/(7 yr) = 10"* '^^ Hz which roughly corresponds to 
the optimal sensitivity of the data of Refs. |2^. To get a first idea of the situation, let us start by considering the 
fiducial value agd = 10"'' corresponding to G/i ~ a/50 = 2 x 10~^, i.e. the traditionally considered type of string 
tensions (which can naturally come from Grand Unified Theories (GUT) and which is most relevant for large scale 
structure formation). We can compare three different GW amplitudes of relevance for the pulsar experiment: (i) the 
amplitude of individual bursts having a recurrence rate f^^'' — l/(7yr) (instead of the 1/yr recurrence rate considered 
above for LIGO or LISA), (ii) the rms amplitude of our confusion background (6.17), and (iii) the rms amplitude of 
the usually discussed full integral (6.15) which includes both rare bursts and overlapping ones. We find 



/i™;P(afid,/r) = 0.503 x 10-", 
/iconfusio„(afid,/r''') = LOl X 10-" , 



^r;'(«fid,/r) = 2.30x 10 



-13 



(6.18a) 

(6.18b) 
(6.18c) 



In computing the integrals (p. 150 
space density of loops, ng, Eq. (157 



and (6.17) we have here used (for better accuracy) an improved estimate of the 



Indeed, numerical simulations indicate that, beyond the scaling ni ~ a^" t^" one 
must add an extra factor related to the parameter characterizing the density of long strings. This factor is different 
in the radiation era and in the matter era. In the notation of |l|] this extra factor is ~ 0.4 ~ 10 in the radiation 
era, and is ~ 0.12 Cm 1 during the matter era. In other words, a better estimate of ne is obtained by multiplying 
the estimate (5T) by the function 



C{z) = l + 9z/(z + Zeq) 



(6.19) 



which interpolates between 1 in the matter era and 10 in the radiation era. 

Note that, in terms of the contribution (per frequency octave) of GW's to the present energy density, pgw(/) 
{2tt f h{f))'^ / (WttG) ~ (tt/AG) p h? , or, better, of their fractional contribution to the closure density. 



Pc ^ 



(6.20) 



the results (3.18) yield 



confusion 



(afid,/D = 2.99x 10- 



(6.21) 



f^g™(«fid, /r) - 1-57 X 10-6 . (6.22) 

The usually considered fig^ is 5.25 times larger than the physically meaningful confusion noise! This large discrepancy 
comes from the fact that fig^"^' includes the time-average contribution of rare, intense bursts, which are in general, 
not relevant for a pulsar experiment (if they are so rare that they do not occur during the actual duration of the 
experiment) . 

From this consideration, it would seem to follow that the usual way to use pulsar data to set limits on G fi (i.e. the 
comparison between the theoretically predicted ng™j^'j5ig(a, /) and the observational constraint on a Gaussian W^{f)) 
is seriously affected by the present work. It would also seem that the correct way to set limits on Gfi fro m pulsar data 



consists simply in replacing Og^"^^j^g by our new, significantly smaller, fig™^"*""". Then the value (6.21) suggests that 



even cusp- (rather than kink-) dominated backgrounds, with c = 1 (as the one considered in the equations above) 
generated by string tensions G fj. oi order of 10"^ might be compatible with pulsar constraints at the Sl°^*'<10"'' 
level. [In view of the crudeness of our estimates we shall not try here to give any precise limit on G /U from fi' 



obs 



However, Eqs. (p. 18) show that the situation is actually somewhat more complex than that. Indeed, Eq. (6.18a) 



shows that the (observationally relevant) 1/(7 yr) bursts have an amplitude comparable to the full confusion noise 



(which sums many overlapping, small signals). [This comes from the fact that the confusion integral (6.17) is dominated 
by its lower limit.] Therefore a significant part of the difference between /iconfusion and /lusuai comes from not very 
intense, but not very rare bursts. [Note, however, that the dominant contribution in /i^m"'^' comes from the very 



intense, very rare bursts, with recurrence time 3> 7yr.] In other words, Eqs. (3. IS) show that, within the frequency 
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FIG. 3. Usual rms noise (upper short-dashed curve), confusion noise (lower short-dashed curve) and burst GW amplitude 
(solid line) emitted by cosmic string cusps, in th e freq uency band relevant for pulsar timing observations (/c — N — l/(7yr)). 
Here, we assumed c = 1 and included the factor ([6. J9|) in the spatial density of loops. 



bandwidth /f"'"' relevant for pulsar timing, the GW signal h{f) is a complicated superposition of a nearly Gaussian 
noise (of variance ^confusion) ^^"^ ^ small number of occasional random bursts, occurring on the {f^^'^)^^ time scale 
and of amplitude comparable to /iconfusion(/c)- In addition, there might also occur (on longer time scales) some larger 
bursts. This situation shows that one needs to reanalyze from scratch the pulsar limits on G by dealing explicitly 
with the statistical properties of such a complicated mix of signals, i.e. by tackling seriously the strongly non Gaussian 
nature (involving an important quasi-Poissonian component) of /i(/) within the pulsar timing bandwidth. Until such 
an analysis, using our new results on the nature of the string GW background, is performed one cannot draw secure 
limits on G II from pulsar observations. We expect, however, that the result of such an analysis will be, to a good 
approximation, equivalent to replacing Jlg^"''^ by our new $7^™^"'''°" (which is about five times smaller than fig^"^' 
when G/i = 2 x 10~^, and about four times smaller when G/x = 10~®). In particular, we expect that our results make 
a GUT-like value G/i ~ 10^^ now compatible (even with many cuspy loops) with present pulsar data. 

Leaving to future work such an analysis, we content ourselves by comparing in Fig. ^ the variation with a of the 
cusp burst signal (for c = 1 and N = 1/(7 yr)), the confusion GW amphtude /iconfusion(/c *''^) ifF^'^ = l/(7yr)), and 
the usually considered total rms amplitude /ir,nT'(/c"''^)- Contrary to what happened in Fig. |l| and Fig. ||, we see 
now that the burst signal and the confusion signal are of comparable orders of magnitude in a wide range of values 
of a — 50 G/i. Note also that /i",™^' is a significant overestimate of /iconfusion when Q!>5 x 10~^, i.e. G/i>10~'', 
which includes the GUT case which has been traditionally of most interest for cosmic string research. In view of the 
subtlety we just mentionned concerning the data analysis of pulsar experiments we did not indic ate in Fig. ^ a precise 
"one sigma" level for the sensitivity of the pulsar experiment. [A rough guess, using Eq. ( 6.20| ) with ^ 10~^ is 
hP""' - 0.5 X 10~i3.] 

Switching from a defensive attitude (pulsar limits on G/i) to an optimistic one (detection of GW by pulsar exper- 
iments), and forgetting for a moment the subtlety of the fact that /i^urst ^ ^confusion ^ assuming that /jconfusion 
gives a good first estimate of the GW amplitude to be compared to the timing precision of pulsar experiments), it is 
striking to notice on Fig. || that /jconfusion very flat function of a, so that a modest improvement in the sensitivity 
of pulsar experiments (due either to a longer time span or to the discovery of an intrinsically more stable pulsar) might 
allow one to detect the confusion noise coming from a string network down to a ~ 10~^ (i.e. G ^ 10~^^). Evidently 
one should keep in mind that Fig. ^ is drawn for an average cusp number c = 1. Assuming a smaller value of c, or 
even c = and considering only the smaller kink signals, will make it much more difficult for pulsar experiments to 
probe the existence of cosmic strings. 



VII. CONCLUSIONS 



We have studied in detail the amplitude, frequency spectrum, waveform and rate of occurrence of the high-frequency 
gravitational wave (GW) bursts emitted at cusps and ki nks o f a cosmological network of oscillating loops. Our main 
tool in studying the waveform has been the factorization (2.23) of the Fourier transform of the emitted GW amplitude. 
This factorization allowed us to conveniently extract from the waveform its physically meaningful, gauge-invariant 
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content. [This is why our waveforms differ from previous results which did not notice the gauge nature of 

the leading terms.] In the time domain these waveforms correspond to power-law "spikes" (x \t — tc\^ {(3 = ^ for 
cusps; and /3 = | for kinks) of linearly polarized GW's, with some smoothing of the center of the spike on time scales 
\t ~ tc\ ^ 0^ Tg, where = i/2 is the loop oscillation period and where 6 is the misalignement between the center of 
the beam and the direction of emission. 

We estimated the rate of occurrence and the distribution in amplitude of the GW bursts emitted at cusps and kinks 
by using a simple model for the cosmic string network. When comparing our results with observations, one should 
keep in mind the simplifying assumptions involved in our model: (i) All loops born at time t were assumed to have 
length £ at with a ^ T G and T ^ 50. It is possible, however, that the loops have a broad length distribution 
n{£,t) and that the parameter a characterizing the typical loop length be in the range T G ^ < a < 10~^. (ii) We 
have also assumed that the loops are characterized by a single length scale, with no wiggliness on smaller scales. 
Short-wavelength wiggles on scales ^ F G ^ i are damped by gravitational back-reaction but some residual wiggliness 
may survive, thereby modifying the amplitude and the angular distributio n of the GW bursts from cusps and kinks, 
(iii) In many of our estimates we assumed the simple, uniform estimate (5J) for the space density of loops. This 
estimate is probably accurate in the matter era but is expected to be too small by a factor 10 in the radiation era 
In Section Vl where the contribution to the confusion background of the radiation era was crucial we corrected 
the estimate (5.1) by including the redshift-dependent factor (6.1£). (iv) Finally, we disregarded the possibility of 
a nonzero cosmological constant which would introduce some quantitative changes in our estimates. As a general 
comment, let us recall that, though we tried to keep the important "27r factors", our estimates have systematically 
neglected factors of order unity. 

In our view the most important astrophysical results of the present investigation are the following: (1) the clear 
recognition of the strongly non Gaussian nature of the string-generated GW background. The slow decrease with 
frequency of the GW burst signals means that there are occasional sharp bursts that stand above the "confusion" 
GW noise made of the superposition of the overlapping bursts. (2) GW bursts from cusps might be detectable by 
the planned GW detectors LIGO /VIRGO and LISA for a wide range of string tensions even if the average number of 
cusps per string oscillation is only 10%. In spite of the argument pT[ | that string reconnection can inhibit cusps (which 
are generic for smooth loops fl^]), we find it plausible that 10% or at least a few % of the loops in the network will 
feature cusps. In view of the crucial importance of the average number of cusps for detection by LIGO we recommend 
that new simulations be performed to determine this quantity. (3) Even if the number of cusps turns out to be very 
small, our estimate of the GW amplitude emitted by the ubiquitous kinks show that the space borne GW detector 
LISA has the potential of detecting GW bursts from kinks in a wide range of string tensions. (4) Finally, we show 
the need of a reanalysis of th e constraints on G fi derived from pulsar timing data. Indeed, for such low frequencies 
the usual estimate, Eq. ( 6.15| ), of the GW stochastic background (which neglects its non Gaussianity) seems to be 
quite inadequate because i t ave rages on very rare, intense bursts. We have introduced a new, more relevant quantity, 
the confusion noise, Eq. ( |6.17 ), which averages only over the overlapping bursts. In the first approximation, we 
expect that the usually derived pulsar timing data limit on a Gaussian stochastic back groun d (often expressed as a 
limit f^^'psR^^lO^'') will entail essentially the same limit on the confusion part, Eq.( |6.17 ), of the GW stochastic 
background, i.e. we expect that the real pulsar limit on G/i will be the weaker constraint f7g™^"*''°'^(G/i) < f^^^psn"- 
Our rather crude approximations do not allow us to transform this relaxed limit in a precise limit on G/i. However, we 
expect that our results render a GUT-like value Gfi ^ 10~^ compatible with pulsar data, even if c 1 (and probably 
easily compatible with G/i ~ 10~^ if c is 10% or less). However, we emphasized that there are still occasional bursts 
that complicate the analysis and call for an improved treatment. Until such a careful analysis is done, together with 
a precise estimate of the number of cusps in a string network, one cannot use pulsar data to set precise limits on G fi. 
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